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Vapor deposition on polycrystalline films can lead to extremely high levels of compressive stress, 
exceeding even the yield strength of the films. A significant part of this stress has a reversible nature: 
it disappears when the deposition is stopped and re-emerges upon resumption. Although the debate 
on the underlying mechanism still continues, insertion of atoms into grain boundaries seems to be 
the most likely one. However, the required driving force has not been identified. To address the 
problem we analyze, here, the entire film system using thermodynamic arguments. We find that 
the observed, tremendous stress levels can be explained by the flux induced entropic effects in the 
extremely dilute adatom gas on the surface. Our analysis justifies any adatom incorporation model, 
as it delivers the underlying thermodynamic driving force. Counterintuitively, we also show that the 
stress levels ‘decrease’, if the barrier(s) for adatoms to reach the grain boundaries are ‘decreased’! 


I. INTRODUCTION 


During the growth of a polycrystalline film on a sub¬ 
strate, the film usually develops a significant amount of 
internal stress. If the film temperature is high enough to 
reach Volmer-Weber (VW) type growth conditions [l|, Q , 
the film stress during deposition follows a compressive- 
tensile-compressive (CTC) evolution, as is indicated with 
stages I, II, and III in Fig. [TJ During stage I, the nucle- 
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FIG. 1: Stress evolution during VW-type film growth 

for copper deposition with a flux of 0.1 A s~^ onto silicon 
oxide at room temperature. It consists of three main stages: 
nucleation (I), coalescence (II), and thickening (III). The de¬ 
position was interrupted 3 times for 350 minutes. The re¬ 
versible stress jumps are indicated by Aai (reproduced from 

i). 


ated islands devel^ a compressive stress due to surface 
tension effects [3,0 • Stage II occurs during film closure 
when the 3D islands coalesce and form grain boundaries 
(GBs). At this stage, the film free energy can be further 
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lowered by ‘GB zipping’, which in turn delivers tensile 
stress d, |3]- Without sufficient mobility (low temper¬ 
ature or high deposition flux) the film remains tensile 
upon further growth. In contrast, a maximum tensile 
stress develops for VW-type growth of high mobility ma¬ 
terials, which occurrence coincides approximately with 
the moment the film closes. From this moment on, the 
stress turns once again towards compressive values (stage 
III) Q. Surprisingly, a significant part of the compres¬ 
sive stress has a reversible nature: upon interrupting the 
deposition flux (see Fig. [T]), the film stress jumps to less 
compressive values and the original compressive stress 
state before interruption is almost fully restored when 
the flux is switched on again. These stress jumps can be 
as large as ~150 MPa and the time constant of the stress 
variation upon resuming the deposition is in the order of 
20 seconds ii! 

In the last 20 years several mechanisms have been 
proposed aiming to explain the observed effects: 1) 
pre-coalescence surface tension continuation combined 
with ongoing grain growth 2) surface roughness de¬ 
velopment during deposition combined with step-step- 
interactions , 3) adatom insertion into GBs pTV o^ , 4) 
interaction of adatoms with surface and each other [9|, 5) 
‘inside bundling - outside grooving’ of GBs 6) depth 
changes in the GB grooves [l^, ... Whereas several of 
these mechanisms rely on kinetically limited processes, 
the GB adatom insertion model suggests that the com¬ 
pressive stress is generated via adatoms that are forced 
into the GBs by the enhanced chemical potential (CP) 
of the surface that is set up by the deposition flux. By 
switching off this flux, the CP should drop, which should 
lead to an outflow of the excess atoms from the GBs and, 
thereby, to a relaxation of the compressive stress [HI- Re¬ 
cent experiments confirmed that GBs are prerequisite for 
the existence of the reversible stress jumps [l^ . However, 
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more questions arise, as the time constant of the stress 
relaxation upon interruption seems to be temperature in¬ 
dependent [13 ■ On the other hand, surface stress effects 
@ are expected to be too low in magnitude [l^ to ex¬ 
plain the reversible stress jumps. While the discussion on 
the mechanism(s) still continues, at a more fundamental 
level, the underlying driving force behind the effect has 
never been addressed. In this paper we derive the mag¬ 
nitudes and the changes of the CP on the surface next to 
the position of the GBs and show that this indeed forms 
the driving force for any adatom insertion model. 

From a thermodynamic point of view, the most fun¬ 
damental question has never been addressed, probably 
due to conceptual difficulties in calculating the CP of 
the surface during the growth; ‘how can a flux (change) 
as low as ~0.1 monolayer per second [ML s“^] lead to 
stress jumps as high as ~150 MPa?’ Our study focuses 
exactly on this question and we show not only that these 
low fluxes can generate such huge stresses, but also that 
the driving force for the stress jumps is decreased, if it 
is easier for the adatoms to diffuse towards and into the 
GBs. 


II. RESULTS 


A. Basic Thermodynamic Description 


To derive our model, it is important to realize that the 
film is under growth conditions and therefore naturally 
not in equilibrium. However, as long as the growth con¬ 
ditions do not change, it can be treated in steady-state, 
like the famous ‘ Growth-Wulff construction’ [I 3 . The 
enhanced surface CP with respect to equilibrium sets up 
an adatom current to steps, which finally leads to the 
film growth. This also means that the CP on the surface 
varies locally and that positions connected to each other 
will try to balance their difference. If atom transport is 
sufficiently active on the time scale of consideration, one 
can approximate adjacent positions to be in equilibrium. 
Therefore, for constant ‘small’ deposition fluxes and the 
absence of kinetic limitations, thermodynamic equilib¬ 
rium can be assumed between the positions on the sur¬ 
face immediately next to grain boundaries (s/GB), the 
grain boundaries (GB), and the grain interiors (g). This 
assumption is further underpinned by the small number 
of total additional atom that have to be incorporated 
in the GBs. For the surface we solve rate equations to 
determine the CP immediately next to the GBs and we 
further treat this position, the GBs, and the grain in¬ 
terior to be in equilibrium. Thermodynamic equilibrium 
certainly does not hold for the transition between the flux 
‘on’ and ‘off’ state, but is justified a few tens of seconds 
after the flux change (see above). Therefore, at constant 
or zero flux, a change of the CP of the surface next to 
the GBs, will finally change the CP of the GBs with the 
same amount, which in turn will change the CP of the 


grains: 


^Ms/GB ~ A/tgb ~ A/ig => A/Ts/gb ~ A^g (1) 

This core equation enables us to bypass the determina¬ 
tion of the CP of the GBs, as well as the absolute CP 
values on the surface and within the grains. 


B. Chemical Potential of the Surface 


The free energy of a surface depends on the formation 
and interaction energies of a myriad of surface features 
such as terraces, steps, kinks, step adatoms, adatoms,... 
As the surface morphology evolves during deposition, the 
population of these features changes accordingly. How¬ 
ever, it is known that the reversible, compressive stress 
can develop within seconds after starting the flux even 
with rates as low as 0.1 ML s“^. Obviously, the popu¬ 
lation of point-like features, like adatoms, step-adatoms, 
and kinks, can (and will) change abruptly upon the ar¬ 
rival of flux on the surface, but extended surface features, 
like terraces and step edges, do not change significantly 
within such short time scales, as they consist of a large 
number of atoms [ 2 ^. For example, the surface roughness 
is directly linked to changes in the appearance, distribu¬ 
tion, and amount of steps and terraces. The CTC be¬ 
havior is usually observed under step-flow growth mode 
conditions, where changes in the surface roughness are 
known to happen very slowly . This means that 

the gradual increase in surface roughness (and hence the 
extended surface features) during deposition will only 
have a long time effect on the surface CP via the Gibbs- 
Thomson relation. Since the reversible stress jumps oc¬ 
cur in a matter of seconds, we safely can ignore these long 
term changes in our analysis. Moreover, in contrast to 
the adatoms that live in a 2D-space on the terraces, the 
step-adatoms and step-kinks are confined to the ID space 
on the step edges. This causes the rate, at which the 
step-adatoms and kinks meet and annihilate each other, 
to be significantly higher than the adatoms on the ter¬ 
races. As a result, the increase in adatom population 
on the terraces, upon starting the deposition, is orders 
of magnitude higher than that of kinks or step-adatoms 
[2l[. The conclusion is that the surface GP change that 
is responsible for the almost instantaneous stress jumps, 
is mainly dominated by a change in the adatom density. 
Gonsequently, we can neglect all other contributions, as 
they would lead only to higher order correction terms in 
determining the surface GP variations: 


A/Xg — A^adatom T 0(A^step adatom j A^kinki A/igtep, -.) 

. ^‘5'adatom 


A 


^Cadatom ^ '^t/adatom int. 


ON 


dN 


-T- 


dN 


( 2 ) 


The first term in the brackets, dUg^dutom/dN, accounts 
for the surface temperature dependent change in aver¬ 
age energy (potential and kinetic) of individual adatoms. 
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Given by the radiation of the evaporator and the kinetic 
energy of the arriving atoms, the increase in film temper¬ 
ature is less than 10 K for Cu, Ag, and An [1^, which 
corresponds to ~2.6 meV per film atom according to the 
classical Dulong-Petit limit of the heat capacity in solids. 
As these materials all have an excellent thermal conduc¬ 
tivity, the surface and the bulk temperatures are virtually 
identical. Since we finally have to compare only the ‘CP 
variation’ of the surface and the grains, we can safely ne¬ 
glect this term, as we would have to add the same value 
to both sides of Eq. [T] 

The second term, 9C/adatom int./^A^, corresponds to 
the interaction energy between the individual adatoms 
given by (combinations of) Van-der-Waals, electrostatic 
(dipole), elastic, and electronic (substrate mediated) ef¬ 
fects . We safely can ignore this term, as STM exper¬ 
iments at ^15 K have shown that the absolute value of 
the interaction energy drops below 0.1 meV for two Cu 
adatoms separated more than 60 A on a Cu(lll) surface 
[ 2 ^. This is equivalent to an adatom density (fractional 
coverage) of < 6.0 x 10“^ ML, and as it will be shown 
in the following, we never reach such densities during the 
deposition. 


3 <inter-adatom dist.> (nm) 
1000 100 10 1 



Adatom density, 9 (ML) 



FIG. 2: Surface CP on terraces and next to GBs: (a) 

entropic component of the CP on Cu(lll), estimated by the 
‘adatom lattice gas’ and the ‘2D adatom gas’ model, (b) 
Adatom density profile, 9„, on a terrace; the corresponding 
energy landscape is indicated at the top. (c) Typical step 
configuration in the vicinity of a GB caused by the ‘Zeno 
effect’ (top). As AAes vanishes for terrace widths with w < 6, 
the combination of these terraces can be approximated with 
one effective terrace (bottom). This leads to a ‘lower’ adatom 
density near the GB and, therefore, to ‘lower’ stresses. 


The third term, —T x dSadatom/dN, involves the en¬ 
tropic effects of the 2D adatom gas. In general, depend¬ 
ing on the adatom mobility, adatoms can be assumed to 
be confined on discreet lattice sites (‘adatom lattice gas’) 
or to be delocalized behaving as a 2D Van-der-Waals sur¬ 
face gas (‘2D adatom gas’) (see Supplementary Note 
1 ). As these two models naturally set the lower and up¬ 
per limits for the adatom gas entropy, we calculated the 
boundary values of the CP for copper in Fig. [2^. Al¬ 
though the absolute values differ more than 0.3 eV, both 
models show a linear behavior in this logarithmic plot be¬ 
low 0.01 ML such that the following approximation holds 


for the CP variations: 


A 



kBT\n{02/ei) 


( 3 ) 


C. Adatom Density on Terraces 

To calculate the adatom densities during deposition 
and interruption. Fig. [21 d shows a simplified model of the 
film surface, in which we define the position of the first 
lattice row next to the ascending step edge as the origin 
of a terrace with width w in lattice units. By solving the 
differential equation for mass conservation, the adatom 
density at site n on the terrace, can be derived as a 
function of deposition flux F (see Supplementary Note 
2 ): 

n /, , Fw{an + l){sw + 2) Fn^ 

2,xd{asw + a + s) 2,,^ 

, where 6»eq = exp (-Fform/fcBT), L'd = 
L'oexp (-Adiff/ZcBT), a = exp (-AFatt/^BT) and 
s = SQexp{-AEEs/kBT), in which uq, Edis, Ffomi, 
AFatt, So, AEes are the diffusion rate prefactor, diffu¬ 
sion barrier, adatom formation energy from a kink site 
of the step, the attachment barrier, correction prefactor 
for hop over the step, and the Ehrlich-Schwoebel barrier 
(see Fig. Hi). Note that at zero deposition flux, the 
adatom density at each position n of the terrace is equal 
to the equilibrium density 0eq- For constant deposition 
with AFes ^ AFatt, s <S a, the adatom density is 
highest close to the end of the terrace (see Fig. m), 
whereas the maximum shifts to the middle of the terrace 
for low values of AEes (see Fig. [S}:). Note that the 
maximum of the adatom density is only exactly at the 
end of the terrace for Fes = 00. 

Combining Eqs. [2] to IH one can calculate A fig as a 
function of deposition flux for any site n on the terrace. 


D. Chemical Potential of the Grains 


Considering an in-plane isotropic biaxial film (dx = Cy 
and dz = 0), it can be proven for the right hand side of 
Eq. [1] that the CP within the grains is proportional to 
its total internal stress level dg, (see Supplementary Note 
3): 


A^g — 


.+A<tL' 


— Dx^^dx '^^dy 


, where D is the atomic volume. 


2 

—flAdg 

( 5 ) 


E. Derivation of the Stress Jumps 

Motivated by the fact that stress emergence has to be 
GB related [l^, it is crucial to derive the stress jumps 
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using the A/^tg immediately next to the GBs! Combining 
the above equations at position n = w, the predicted 
upper limit for the reversible stress jumps based on pure 
thermodynamics is given by 


(I _|_ Fw{aw + 2) \ 

2 n \ 2i/d (a s w + s + a) 0eq J 


( 6 ) 


It has been shown for gold at room temperature that 
step-flow growth is the underlying atomic process for 
polycrystalline film growth in the VW regime [ 1 ^. On 
the basis of the so-called ‘Zeno effect’, terraces closer to 
the GB get progressively decreased in their width during 
depostion, which leads to an enhanced surface curvature 
at the GB vicinity, as sketched in Fig. [2t-top [13, H^l- 
This results in a deviation from the macroscopic equi¬ 
librium surface of an annealed polycrystalline film |28l| . 
Although not mentioned in [ 13 ], this deviation changes 
only slightly back over one hour upon stopping the depo¬ 
sition while keeping the film at room temperature. The 
same effect, which is due to the existence of a signifi¬ 
cant Ehrlich-Schwoebel barrier, has also been observed 
on Gu(lll) [13|- For our study, we, therefore, safely 
omit Gibbs-Thompson correction terms associated with 
macroscopic surface curvature variations. 

The red dashed lines in Fig. [3] show the predicted 
stress jumps derived via Eq. [5] for surface terrace widths 
between 1 and 500 atomic spacings and the existence of 
an Ehrlich-Schwoebel barrier. It is striking that the ob¬ 
tained stress values exceed even the experimentally ones 
(crosses) and that we receive a rather good agreement 
already for a terrace width with w = 1 atomic spacing. 
To get a feeling for the numbers, the top horizontal axis 
shows the adatom density at the end point {w = n) for 
a terrace with a width of rc = 500 (74 nm). Note that 
the adatom density is less than 10“^ ML even for fluxes 
as high as 100 ML s“^. This validates our dilute adatom 
gas assumption while calculating the surface GP and jus¬ 
tifies our steady-state thermodynamic approach. Note 
also that it is surprising that a dilute adatom gas of less 
than 10“^ ML has the potential to induce '^1 GPa stress 
in the film, which is more than both the yield and the 
ultimate strength of copper! Please note that this stress 
can be realized in the bulk only, if there exists a kineti- 
cally not limited atomic mechanism that transfers the GP 
variation of the surface to the grain interior. The curves 
are calculated for Gu(lll) and we have used T = 298 
K, = i(361.49 pm)3, i/q = Hz^diff = 0.040eV 
m So = 15 [ 13 , AEes = 0.224 eV [H, AE^^t = 0 eV 

(AEatt ~ 0 for most metals at room temperature), and 
Eform = 0.714 eV [111, for our calculations, which implies 
an adatom equilibrium density of 0eq = 8.6 x 10“^^ ML 
at zero deposition flux. 

It is known that stress relaxation mechanisms are ac¬ 
tive both during the growth and after stopping the de¬ 
position, which reduce the absolute intrinsic film stress 
[3 a [ 3 ^ . Indeed, the experimental stress values (black 
crosses) are in general lower than the red dashed stress 
lines. However, if one considers that the time constant 


Adatom density (ML) 

at the end (n=500) of a terrace with a width of w=500 



10 10 10 
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FIG. 3: Reversible stress jumps for a (111) textured 
copper film as a function of flux: The dashed red 
curves are calculated for typical terraces (!<«;< 500 and 
AEes = 0.224 eV), whereas the solid blue curves describe 
terraces in the vicinity of GBs {wes < 20 and AEes ~ 0 eV 
) where funneling takes place. The top horizontal axis shows 
the adatom density for a typical surface terrace {w = 500 and 
AEes = 0.224 eV). Crosses show experimentally reported 
literature values: 1,2,3 [3, 4 [^, 5,6 |3], 7 [§3, and 8 [^ . 
Note that the stresses are ‘lower’, if it is ‘easier’ for atoms to 
diffuse to and into GBs. 


of the stress relaxation processes are distinctively larger 
than the time constant of the reversible jumps [^, the 
discrepancy between the observed and the predicted val¬ 
ues for typical surface terraces is too large to be explained 
by the stress relaxation effects alone. 


F. Stress Jumps considering Funneling 


As our derivation of the stress jumps exceeds in gen¬ 
eral the experimental values, we turn our attention to 
experiments on both Gu(lll) and Ag(lll) at room tem¬ 
perature: these experiments revealed that the Ehrlich- 
Schwoebel barrier of the lower step vanishes,if the dis¬ 
tance between two neighboring steps becomes less than 
6 atomic spacings [^ [^ ! This effect opens a fast mass 
transfer channel for terraces with w < 5 called ‘funnel¬ 
ing’. The fast mass transport region can extend up to 
21 atomic sites (over 5 steps) away from the GB. We 
account for this by setting the critical terrace width to 
w = 5 and by assuming that the Zeno effect reduces the 
width of the subsequent terraces by one atomic spacing. 
With this approximations, the last 6 terraces next to a 
GB can be treated as one single terrace with an ‘effective’ 
width of Weff — 20, with AEes = 0 (s = 1), see Fig. O: 
bottom. However, as intermediate step edges given by 
the 5 steps can potentially act as adatom sinks, we eval¬ 
uate the stress jumps for effective terrace widths, WeS, 
ranging 1 to 20 atomic spacings (see blue lines in Fig. 
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m- As our calculations define upper limits for the stress 
jumps, a tunneling width of 5 spacings represents the 
best fit. This means that is it enough, if only the last 2 
steps (and not 5) before the GB show funneling. The fact 
that our result with the inclusion of funneling delivers a 
rather good Ht with the experimental values is a strong 
indication for the validity of the GB insertion model, es¬ 
pecially, as we calculate the CP on the surface exactly 
next to the GBs. Note, however, that we do not address 
the exact atomic pathway (e.g. diffusion, exchange,...), 
as we evaluate only the thermodynamics. 


III. DISCUSSION 

Although we ‘lowered the barriers’ for the atoms to dif¬ 
fuse towards and into the GB, we receive ‘lower stresses’ 
than without funneling. This seemingly counterintuitive 
behavior demonstrates that we are not addressing a cer¬ 
tain atomic diffusion/incorporation model, but calculate 
the thermodynamic driving force on the basis of the CP. 
Lowering the barriers for atoms to diffuse towards and 
into the GB, decreases the adatom density near the GB 
and results in a reduction of the driving force for the 
atoms to diffuse into the GB. Evidently, the funneling 
curves predict the correct order of magnitude for the 
stress levels. By setting (s = a = 1) in Eq. SI one gets 
an estimate of the adatom density at the GB vicinity 
n = WeS- for a deposition rate of 1 ML the adatom 
density is predicted to be lower than 10“^° ML! 

Finally we address another hypothetical mechanism in 
combination with the equilibrium situation between the 
surface and the bulk. For materials with a low enough 
Ehrlich-Schwoebel barrier and funneling terraces, the CP 
distribution shows a maximum around the middle of the 
terrace, see Fig. 2c. In addition, the CP at the end posi¬ 
tion of, e.g., a large terrace in the middle of the grain 
(with Ehrlich-Schwoebel barrier) is significantly larger 
than the CPs of funneling terraces that connect to the 
grain boundaries. Pure thermodynamic considerations 
imply that also these maxima tend to establish equilib¬ 
rium with the interior of the grain underneath. As the 
grain is comparable to a single crystal, dislocation nu- 
cleation would be an imaginable pathway to balance the 
CP differences. However, on epitaxial films and single 
crystals the reversible stress jumps are not observed [1^. 
The reason for this is a high nucleation barrier: a critical 
stress of 1.3 GPa has been determined to nucleate dislo¬ 
cations in Cu at 300 K [d^. Such high absolute stress 
values are neither observed experimentally, nor does our 
model predict an equivalent rise of the surface CP (except 
for large terraces in combination with high deposition 
fluxes). Dislocation nucleation is, therefore, kinetically 
limited! The overall picture is that during the growth 
all points on the surface are in parallel trying to balance 
their CPs with adjacent positions as well as with the 
grain underneath. However, in which way and by which 
rate the CP of the grains will change clearly depends on 


the rates of the pathways between all subsystems: sur¬ 
face, GB, and bulk. With a significant dislocation nu¬ 
cleation barrier and active GB diffusion as well as atom 
incorporation, the whole system quickly evolves towards 
equilibrium via atom insertion in GBs. As a result the 
CP difference between the grain and the surface CP max¬ 
ima are reduced, which effectively lowers the driving force 
for dislocation nucleation making this latter process even 
less favorable. 

If one intends to compare our results with experiments, 
it is important to realize that we determine only the 
pure reversible equilibrium jumps. For a proper com¬ 
parison the experiments should have no kinetic limi¬ 
tation of atoms going in/out of the grain boundaries, 
should be performed long enough such that equilibrium 
has reached (all GBs show the equilibrium density of 
additional atoms), and no stress relaxation mechanisms 
should occur. In this limit, we expect the stress jumps 
to be GB density independent. Kinetic limitations would 
immediately result in a GB density dependence, as equi¬ 
librium will not be reached and the rate towards equilib¬ 
rium scales with the number of the pathways and hence 
the GB density. 

The deposition flux and temperature dependence is more 
complex, as the growth mode (layer-by-layer, step flow, 
3D/rough growth) that determines the size of the ter¬ 
race next to the GBs, also changes with deposition rate 
and mobility. If one, e.g., lowers the rate for a film that 
growths in 3D growth mode, one might enter step flow 
growth conditions in which effective larger terraces (with 
higher adatom density) might communicate with the GBs 
such that the stress jumps are even higher instead of 
lower. 

Our analysis shows that entropic effects in the ex¬ 
tremely dilute adatom gas on the surface of a polycrys¬ 
talline film during vapor deposition are strong enough 
to cause plastic deformation in the film. The predicted 
film stresses are even higher than the observed ones. 
If we ‘lower the barriers’ for atoms to diffuse towards 
and into the GB by funneling, the ‘stresses decrease’ 
and the predicted values perfectly match the experimen¬ 
tal ones. With this we deliver the, until know missing 
thermodynamic driving force for any GB atom insertion 
model. Further experimental research, similar to [T5l.[^. 
is needed to clarify the exact atomistic mechanisms and 
pathways behind this effect. 
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